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Abstract 


The minimum weight design of helicopter 
rotor blades with constraints on multiple 
coupled flap-lag natural frequencies has been 
studied in this paper. A constraint has also 
been imposed on the minimum value of the 
autorotational inertia of the blade to ensure 
sufficient rotary inertia to autorotate in 
case of an engine failure. A stress 
constraint has been used to guard against 
structural failure due to blade centrifugal 
forces. Design variables include blade taper 
ratio, dimensions of the box beam located 
inside the airfoil and magnitudes of the 
nonstructural weights. The program CAMRAD has 
been used for the blade modal analysis and the 
program CONMIN has been used for the 
optimization. In addition, a linear 
approximation involving Taylor series 
expansion has been used to reduce the analysis 
effort. The procedure contains a sensitivity 
analysis which consists of analytical 
derivatives of the objective function, the 
autorotational inertia constraint and the 
stress constraints. A central finite 
difference scheme has been used for the 
derivatives of the frequency constraints. 
Optimum designs have been obtained for both 
rectangular and tapered blades. Using the 
method developed in this paper, it is possible 
to design a rotor blade with reduced weight, 
when compared to a baseline blade, while 
satisfying all the imposed design 
requirements . The paper also discusses the 
effect of adding constraints on higher 
frequencies and stresses on the optimum blade 
weight and the distributions of mass and 
stiffness in the optimum designs. 
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box beam height variation along 
blade span 
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box beam wall thicknesses 
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box beam cross sectional area 
autorotational inertia 
Young' s modulus 
objective function 
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torsional stiffness 

total principal area moments of 
inertia about reference axes 
length of j th segment 

total mass of j th segment 

total number of blade segments 
number of design variables 
blade radius 
total blade weight 

blade weight as a function of design 

variable 0 

box beam weight 

nonstructural blade weight (weight 

of skin, honeycomb, etc. along with 
tuning/lumped weights) 
prescribed autorotational inertia 
design variable increment 
taper ratio in z direction 

i th design variable 

mass density of the j t * 1 segment 

stress in segment 
maximum allowable stress 
blade RPM 


Subscripts and Superscripts 

r root value 

t tip value 

L lower bound 

U upper bound 

A approximate value 
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Introduction 

Computer-based mathematical programming 
methods for optimum design of structures have 
been under rapid development during the last 
two decades. Using mathematical processes, 
engineering design synthesis problems can be 
posed as sequences of analysis problems 
combining engineering models with minimization 
techniques . An extensive amount of work has 
been done in developing such design 
optimization procedures over the past few 
years to bring the state of the art to a high 

level 1-5 . These methods can now be applied to 
optimum design of practical structures such as 

aircraft 1, 5 and helicopters 2-5 . The present 
paper focuses on helicopter rotor blade 
design . 

The helicopter rotor blade design process 
requires a merging of several disciplines, 
including dynamics, aerodynamics, structures, 
and acoustics. Two of the major criteria in 
rotor blade design have been low weight and 
low vibration. For a helicopter in forward 
flight, the nonuniform flow passing through 
the rotor causes oscillating airloads on the 
rotor blades . These loads in turn are 
translated into vibratory shear forces and 
bending moments at the hub. One important 
design technique is to separate the natural 
frequencies of the blade from the harmonics of 
the airloads to avoid resonance. Failure to 
consider frequency placement in the predesign 
stage of the design process could cause a 
significant increase in the final blade weight 
since it generally involves postdesign 
addition of nonstructural masses. To avoid 
such weight penalties it is desirable in the 
design and fabrication of the blade to 
appropriately place the natural frequencies at 
an early stage in the design process. This 
can be done by a proper tailoring of the blade 
mass and/or stiffness distribution. This 
tailoring is not an easy task because of the 
complicated vibration modes of the blade due 

to the presence of several coupling effects®. 
One such coupling is between flap, lag, and 
torsional motions through the pitch angle 
blade twist and offset between the elastic and 
inertia axes. The inclusion of these coupling 
effects makes the design process highly 
complex. In the past, the conventional design 
process was controlled mainly by the 
designer's experience and the use of trial and 
error methods. 

Today, one of the more promising 
approaches to the helicopter rotor design 
process is the application of optimization 
techniques. A considerable amount of work has 
been aimed at optimum designs of vibrating 
structures. For example, minimum weight 
designs with constraints on natural 
frequencies have been addressed in Refs. 7-9 
and the dual problem of maximizing the 
frequencies with a constraint on the total 
weight has been addressed in Ref. 10. 
Frequencies of coupled bending-torsion modes 
caused by an offset between the elastic and 
inertia axes have been addressed in Refs. 9 
and 10. Recently there have been a number of 
applications of optimization techniques to 

rotor blade design 5 ' 6 ' 11-19 . Some of this 
work has been devoted to reducing vibration by 
controlling the vertical hub shears and 

moments 12-17 . In Ref. 13 Taylor described the 


use of modal shaping. The objective of his 
work is to reduce vibration levels by 
modifying 'modal shaping parameters' which are 
functions of blade mass distributions and mode 
shapes. These modal shaping parameters have 
been sometimes interpreted as 'ad hoc' 

optimality criteria 15 ' 17 . In Ref. 14 3ennett 
described a method for reducing the vertical 
shear transferred from the rotor blade to the 
mast by combining conventional helicopter 
engineering analysis with a nonlinear 

programming algorithm. Friedmann* 5 considered 
the problem of minimizing hub shears or hub 
vibratory rolling moments subject to 
aeroelastic and frequency constraints. An 
early attempt at optimum blade design for 
proper placement of natural frequencies with a 
constraint on autorotational inertia was due 

to Peters 16 where he started with a baseline 
blade design and attempted to refine the 
design by trying to find a mass and stiffness 
distribution to give the desired frequencies. 
Reference 17 addressed the optimum design for 
a typical soft in-plane hingeless rotor 
configuration for minimum weight using 
optimality criteria approach. The results in 
Ref. 17 indicate that application of 
optimization techniques leads to benefits in 
rotor blade design not only through 
substantial weight reduction but also a 
considerable reduction in the vibratory hub 
shears and moments at the blade root. In Ref. 
18, Peters addressed a problem of the optimum 
design of a rectangular blade for proper 
placement of frequencies. However, he did not 
use the blade weight as the objective function 
due to a difficulty in finding a feasible 
initial design. Rather, he started his design 
with an objective function involving measures 
of the closeness of frequencies to desirable 
frequencies . 

Currently at the NASA Langley Research 
Center, there is an effort to integrate 
several technical disciplines in rotorcraft 
design. The present work is part of this 
effort and deals with the dynamics aspect of 
design. The problem addressed in this paper 
is an extension of the problem addressed by 
the authors in Ref. 19 where constraints were 
imposed on the first lead-lag dominated mode 
and the first flapping dominated mode along 
with a rotary inertia constraint to assure 
that the blade could autorotate. The 
structural safety of the design was included 
as a first approximation by imposing lower 
bounds on the structural design variables. 
However, the danger of the higher frequencies 
falling in the critical ranges and causing 
resonance remained. The current work involves 
minimum weight designs of helicopter rotor 
blades subject to the following constraints: 
a) upper and lower bounds ('windows') on 
multiple adjacent natural frequencies, b) 
minimum prescribed value on the blade 
autorotational inertia, and c) upper limit on 
the blade centrifugal stress. In Ref. 18 
Peters addressed the necessity of using a 
stress constraint in frequency placement 
optimization but did not include it in the 
optimization formulation. The expression for 
the stress presented here differs from that of 
Ref. 18 and is a more conservative estimate. 
An existing adequate blade which will be 
referred to as the 'reference blade' has been 
selected. In rotor blade design it is 
essential for natural frequencies to be 
separated from values which are certain 
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integer multiples of the rotor speed to 
resonance. These critical values are referred 
to as ' n per rev' where n denotes the total 
number of blades . A modal analysis of the 
reference blade showed that the frequencies of 
interest were away from the n per rev values. 
Hence, it was dacided to define the frequency 
constraints to force the frequencies to be 
close to those of the reference blade. This 
is done by optimally tailoring the blade 
stiffness and mass distributions by the 
procedure developed in this paper. The 

purpose of this paper is to describe the 
formulation and implementation of the 

optimization procedure, present results from 
the procedure and assess the effects of 
additional frequency and stress constraints on 
the optimum designs. 
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the stress due to centrifugal 
is the cross sectional area of 


Optimization Problem Formulation 


The purpose of the optimization procedure 
is to reduce the weight of a blade while 
constraining the natural frequencies to be 
within the 'windows' of the reference blade 
frequencies. The concept of 'windows' has 
been used since the nonlinear programming 
method used in this work cannot handle 
equality constraints. These windows are on 
the frequencies of the first three lead-lag 
dominated modes and the first two flapping 
dominated modes (elastic modes only) . A 
prescribed lower limit on the blade 
autorotational inertia and an upper bound on 
the blade centrifugal stress have also been 
used. Side constraints have been imposed on 
the design variables to avoid impractical 
solutions. The design variables include box 
beam dimensions, taper ratio and magnitudes of 
the nonstructural weights located inside the 
box beam. The optimization process begins 
with an arbitrary set of design variable 
values . 


the i th segment, Mj is the total mass of the 

j'"*’ 1 segment and Cl is the blade RPM. The 
frequencies associated with the first five 
elastic modes of coupled vibration are denoted 
by f^, f 2 , f 2 , f 4 and f 5 , (includes three 

lead-lag and two flapping) . 

The optimization problem can now be 
mathematically posed as follows: 

minimize W ( ) 


where the weight W is given by equation (2) 
and $ denotes the vector of design variables, 
subject to the normalized constraints 


g lc ( < K) - <f k /f ku > - 1-0 

“ 1 - <V f k L > * 0 

g xl (<t>) - 1 - (Al/a) £ 0 

9 11+ k<*> " 1 - < W«V FS > 


k-1,2, . .,5 (5) 
k-1,2, . .,5 (6) 

(7) 

£ 0 k-1, 2, . . , N (8) 


The blade weight, W, has two components 
as follows: 


and side constraints 
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where denotes the box beam weight and W Q 

represents the nonstructural weight of the 
blade which includes the weight of the skin, 
honeycomb, etc., along with the weight of the 
tuning/lumped masses added to the blade. The 
blade is discretized into finite segments and 
the blade weight in discretized form is given 
below: 

N N 

W - £ PjVi + £ W °j <2! 

j-1 j-1 

where N denotes the total number of segments 
and p.. A., Lj and W denote the density, the 

J J J y j 

cross sectional area, the length and the 

nonstructural weight of the segment, 

respectively . 

The autorotational inertia (AI) of the 
blade is calculated as follows 


In equations (5) and (6), f. and f, , denote 

K U 

the upper and lower bound on the k 1 "*" 1 frequency 
f k . In equation (7) a represents the minimum 

prescribed autorotational inertia value. In 
equation (8) O k is the stress in the k th 
segment given by equation (4) , ° max is the 

maximum allowable stress in the blade and tS 
is a factor of safety. In equation (9) <6,. 

denotes the i crl design variable and 0. and 

-u 

0. represent the associated upper and lower 

bounds, respectively. By convention a 
constraint g(0) is satisfied when g(0)£O. 

Analysis 

The modal analysis portion of the program 
70 

CAMRAD which uses a modified Galerkin 
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over three 


2 1 

approach has been used. According to Ref. 
22, this approach is the preferred method for 
computing mode shapes and frequencies of 
structures having large radial variations in 
bending stiffness. Analytical expressions 
have been obtained for the derivatives of the 
objective function, the autorotational inertia 
constraint and the stress constraints. A 
central difference scheme has been used for 
the derivative of the frequency constraints 
(initial attempts using a forward difference 
scheme gave highly inaccurate derivatives) . 



The basic algorithm used is a combination 
of the general-purpose optimization program 

CONMIN and piecewise linear approximations 
for computing the objective function and 
constraints. Since the optimization process 
requires many evaluations of the objective 
function and constraints before an optimum 
design is obtained, the process can be very 
expensive if full analyses are made for each 

function evaluation. However, as Miura 9 
pointed out, the optimization process 
primarily uses analysis results to move in the 
direction of the optimum design; therefore, a 
full analysis needs to be made only 
occasionally during the design process and 
always at the end to check, the final design. 
Thus, various approximation techniques can be 
used during the optimization to reduce costs. 
In the present work, the objective function 
and constraints are approximated using a 
piecewise linear analysis that consists of 
linear Taylor series expansions for the 
objective function and the constraints based 
on the design variable values from CONMIN and 
the sensitivity information from the full 
analysis. Specifically, if the objective 
function F, the constraint g, and their 
respective derivatives are calculated for the 
design variable $ k using an exact analysis, 

their values for an increment in the design 
variable A$ k are as follows: 

NDV 

F = F + (aF/a^jA^ (10) 

k-1 

and 

a ucy 

g - g + T' <ag/a<|> k )A<D k (11) 


where the quantities denoted (~) represent 
approximate values and NDV denotes the number 
of design variables. The assumption of 
linearity is valid over small increments in 
the design variable values and does not 
introduce large errors if the increments are 
small. Since the objective function and the 
constraints are all linearized, the 
optimization problem reduces to essentially a 
sequential linear programming problem. 

A flow chart describing the optimization 
procedure is shown in Fig. 1. The iteration 
scheme is stopped when the objective function 
converges . For the convergence of the 
objective function, a change within a 


convergence tolerance of 0.5 X 10 ^ 
consecutive cycles has been allowed. 


Test Problem 


The reference blade (Refs. 18-19) shown 
in Fig. 2 is articulated and has a rigid hub. 
The blade has a rectangular planform, a 
pretwist and a root spring which allows 
torsional motion. The box beam with unequal 
vertical wall thicknesses is located inside 
the airfoil. As in Ref. 19, it is assumed 
that the box beam contributes to the blade 
stiffness and the contributions of the skin, 
honeycomb, etc. to the blade stiffness are 
neglected. The details for calculating the 
box beam section properties can be found in 
the Appendix of Ref. 19. The properties of 
the box beam located inside the airfoil (Fig. 
2) are as follows: 

h-0.117 ft 
b=0.463 ft 
p»8 .645 slugs/ft 3 
£=2.304 x 10 9 lb/ft 2 

An allowable stress ° max ”l • 93x10 ' lb/ f t 2 and a 

factor of safety FS=3 have been used in the 
analysis. The blade has been discretized into 
ten segments and details of the blade segment 
data are presented in Table 1. The entry 
'min. nonstructural segment weight' in Table 1 
represents the weight of the skin, honeycomb, 
etc. of a segment and 'total nonstructural 
segment weight' represents the weight of the 
skin, honeycomb, etc. along with the 

lumped/tuning weight of that segment. The 
rotor preassigned parameters (the parameters 
that remain fixed during the optimization 
process) are presented in Table 2. 

The frequencies of interest of the 

reference blade are presented in Table 3. The 
first three lead-lag dominated and the first 
two flapping dominated modes are away from the 
critical frequencies (e.g., 3, 4, 5 and 8 per 
rev) and need not be improved further. 
Therefore, the frequency windows for the 

optimum blade are set to be within +1 percent 
of these values (Table 3) . 

Blades with both rectangular and tapered 
planforms have been considered. In case of 
the rectangular blade, the box beam is uniform 
along the blade span. For the tapered blade 
it is assumed, as in Ref. 19, that the box 
beam is tapered (Fig. 3a) and the additional 
design variables are the box beam height at 
the root, h,., and the taper ratio, X^, which 

is defined as the ratio of the box beam height 
at the root to the corresponding value at the 
tip (Fig. 3a) . As in Ref. 19, a linear 
variation of the box beam height, h, in the 
spanwise direction (z direction) has been 
assumed (Fig . 3b ) . 


Results and Discussion 

This section of the paper presents 
results obtained by applying the optimization 
procedure, described previously, to the 
optimum design of both rectangular and tapered 
rotor blades. First, optimum designs are 
described and compared with the reference 
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blade (Refs. 18 and 19). Second, £ftl?PQF aQUALITYrder to satisfy the rotary inertia 


study assessing the effects of including 
constraints on higher frequencies are 
described. Finally, the effects of stress 
constraints have been investigated by 
comparing the results obtained with stress 
constraints to those without stress 
constraints. A summary of the cases studied 
is given in Table 4. Results of these studies 
are presented in Tables 5 and 6 for the 
rectangular blade (30 and 40 design variables) 
and in Table 7 for the tapered blade (42 
design variables). In each table, column 1 
represents the reference blade data; column 2 
gives the corresponding information for the 
optimum design with constraints on the five 
frequencies, autorotational inertia and 
stress (case a. Table 4) ; column 3 gives 
results for the optimum design with 
constraints on the five frequencies and 
autorotational inertia only (no stress 
constraints, case b, Table 4) and column 4 
presents the results (Ref. 19) for the optimum 
design with constraints on the first two 
frequencies (elastic modes only) and 
autorotational inertia (case c, Table 4) . In 
all cases convergence typically has been 
achieved in 8-10 cycles. 

The tables indicate that with the 
constraints on the five frequencies, the 
autorotational inertia and the blade stresses, 
the optimum rectangular blade is 2.67 to 4.74 
percent lighter than the reference blade and 
the optimum tapered blade is 6.21 percent 
lighter than the reference blade. The first 
lead-lag frequency (f^) is at its prescribed 

upper bound after optimization and the 
autorotational inertia constraint is active 
(i.e. exactly satisfied) in all cases. The 
associated design variable distributions are 
presented in Figs. 4-6. Fig. 4a presents the 
optimum versus the reference blade box beam 
horizontal wall thickness ( t x ) distributions 

along the blade span for the rectangular blade 
and Fig. 4b presents the same for the tapered 
blade. In both cases, the optimum blade has a 
larger value of t 1 than the reference blade at 

the blade tip and in case of the tapered blade 
the value of t.^ at the blade root is much 

smaller than the value for the reference 
blade. Figs. 5a and 5b present the optimum 
versus the reference blade box beam vertical 
wall thickness ( t 2 ) distributions along the 

blade span for the rectangular and the tapered 
blade, respectively. The optimization process 
does not produce significant changes in the t 2 

distribution for the rectangular blade (Fig. 
5a) . The changes are more significant for the 
tapered blade (Fig. 5b) where there are larger 
values of t 2 towards the blade tip. The 

larger design variable values towards the 
blade tip are caused by the presence of the 
autorotational inertia constraint which 
encourages the addition of mass at locations 
outboard. Figs. 6a and 6b depict the optimum 
versus the reference blade nonstructural 
segment weight distributions along the blade 
radius. For the rectangular blade (Fig. 6a) 
the optimum blade has lower nonstructural 
weight throughout the blade span. However, 
for the tapered blade (Fig. 6b) the optimum 
blade has larger nonstructural weight towards 
the blade tip than the reference blade. This 
is because the blade is tapered , and has 
reduced structural weight at the blade tip and 


constraint, the nonstructural weight at the 
tip must increase. 

Effect of Constraints on Higher Frequencies 

This section of the paper investigates 
the effect of higher frequency constraints on 
the optimum blade weight and the optimum 
design variable distributions. Therefore, the 
results of the current work which involves 
constraints on five frequencies and rotary 
inertia (case b, Table 4) are compared with 
the results obtained by the authors in Ref. 19 
with constraints on two frequencies and the 
rotary inertia (case c, Table 4) . The results 
of this study are summarized in the last two 
columns of Tables 5 and 6 for rectangular 
blade (30 and 40 design variables, 
respectively) and Table 7 for tapered blade 
(42 design variables) . Table 5 indicates that 
for the rectangular blade with 30 design 
variables, the optimum blade weight increases 
from 89.92 lbs in the two frequency case to 
95.28 lbs in the five frequency case. 
However, the optimum blade with five frequency 
constraints is still 3 percent lighter than 
the reference blade. Tables 6 and 7 indicate 
similar trends for the rectangular blade with 
40 design variables and the tapered blade with 
42 design variables. There is also a change 
in the value of the taper ratio from 1.1 to 

1.5 as shown in Table 7 suggesting that the 
blade taper increases with an increase in the 
number of frequency constraints. The 
optimization process raises the frequency f^ 

(first lead-lag) to its prescribed upper bound 
and the autorotational inertia constraint is 
active in all the cases. 

Figures 7-9 depict the design variable 
distributions (optimum versus reference) for 
the five and two frequency constraint cases. 
Fig. 7a depicts the horizontal box beam wall 
thickness (t L ) distributions along the blade 

span with 30 design variables for the 
rectangular blade. Fig. 7b depicts the same 
distribution for the tapered blade with 42 
design variables. There are significant 
redistributions of the wall thicknesses 
between the five and two frequency constraint 
cases. For example, for the rectangular blade 
(Fig. 7a), in the five frequency constraint 
case (case b) the wall thickness (t^) is 

smaller in magnitude at the blade root than 
the reference blade value but larger than the 
two frequency constraint case (case c) . 
However, at the blade tip the value of t^ in 

the five frequency constraint case is 
significantly smaller than its value in the 
two frequency constraint case, although both 
these values are larger than the reference 
blade value. The situation differs at the 
blade root in the tapered blade case (Fig. 7b) 
where the value of t^ in the five frequency 

constraint case is smaller than the value for 
the two frequency constraint case. Figs. 8a 
and 8b present the box beam vertical wall 
thickness ( 1 2 ) distributions along the blade 

span for rectangular and tapered blades, 
respectively. The figures show that the value 
of t 2 in the five frequency constraint case is 

larger at the blade root than it is in the two 
frequency constraint case whereas the 
tendencies are reversed at the blade tip for 
both the rectangular and tapered blades. 
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Figs . 9a and 9b depict the nonstructural 
segment weight distributions along the blade 
span for the rectangular and the tapered 
blades, respectively. There is a significant 
reduction and change in the nonstructural 
weight distribution between the reference 
blade and the optimum blade in the two 
frequency constraint case than it is in the 
five frequency constraint case. In other 
words, the nonstructural weight distributions 
for the five frequency constraint case is 
closer to that of the reference blade. This 
is because the reference blade was designed 
with a larger number of design requirements on 
frequencies. There are significant 
differences in the optimum design variable 
distribution along the blade span between the 
two and five frequency constraint cases . This 
can be explained as follows. The mass and/or 
stiffness distribution tends to follow the 
pattern of the coupled mode shapes in the 
frequency constrained optimization. In the 
two frequency constraint cases, therefore, the 
mass distributions followed the mode shapes of 
the coupled first lead-lag dominated frequency 
and the first flapping dominated frequency. 
In the five frequency constraint cases, the 
mass distributions followed a different 
pattern as higher coupled frequencies are 
included . 

Effect of Stress Constraints 

The effect of adding centrifugal stress 
constraints to the optimum design with 
frequency and autorotat ional inertia 
constraints has also been investigated. The 
optimum designs with and without constraints 
on the stresses are compared in Tables 5-7. 
Table 5 indicates that for the rectangular 
blade with 30 design variables, the optimum 
blade weight increases with the addition of 
stress constraints. For example in Table 5, 
the blade weight reduction decreases from a 
value of 3.04 percent in case b to a value of 
2.67 percent in case a. The differences in 
weight become more pronounced with an increase 
in the number of design variables (Tables 6 
and 7) . For the tapered blade there is very 
little change in the taper ratio. In all the 
cases studied, the optimization process still 
moves the first lead-lag frequency f^ to it3 

upper bound and the autorotational inertia 
constraint remains critical. 

Some typical results showing the effect 
of stress constraints on optimum versus 
reference blade design variable distributions 
study are presented in Figs. 10-11. Figs. 10a 
and 10b depict the box beam horizontal wall 
thickness (t^) distributions along the blade 

span with and without the stress constraints 
for both rectangular and tapered blades, 
respectively. The presence of stress 
constraints increases the wall thicknesses at 
the blade tip and reduces them inboard for the 
rectangular blade with 30 design variables 
(Fig. 10a) . However, the tendencies are 
reversed in the tapered blade (Fig. 10b) . 
Figs. 11a and lib show the nonstructural 
weight distributions along the blade span for 
the rectangular and tapered blades, 
respectively. For the rectangular blade (Fig. 
11a) the optimization process reduces the 
nonstructural weights at each segment (case b) 
and the inclusion of stress constraints (case 
a) only increases them a little. However for 
the tapered blade (Fig. lib), the stress 
constraints increase the nonstructural segment 


weight at each segment making them higher than 
the reference blade values towards blade 

outboard . 


Concluding Remarks 


In this paper a procedure has been 
described for the minimum weight design of 
helicopter rotor blades with constraints on 
multiple coupled flap-lag natural frequencies, 
autorotational inertia and centrifugal stress. 
The design variables used are the box beam 
cross sectional dimensions, the magnitudes of 
the nonstructural segment weights and the 
blade taper ratio. The program CAMRAD has 
been used to calculate the mode shapes and 
frequencies of the blade and the program 
CONMIN has been used for the optimization. In 
addition, a linear approximation technique 
involving Taylor series expansion has been 
used to reduce analysis time. A sensitivity 
analysis consisting of analytical derivatives 
of the objective function, the autorotational 
inertia constraint and the stress constraints 
and a central finite difference scheme for the 
derivatives of the frequency constraints has 
been performed. Optimum designs have teen 
obtained for blades with both rectangular and 
tapered planforms and compared with an 
existing (reference) blade. Studies have also 
been performed to assess the effects of higher 
frequency constraints and stress constraints 
on the optimum blade designs. 

The following conclusions have been 
drawn from the present study. The 
optimization program CONMIN along with the 
linear approximations based on Taylor series 
expansions has been very efficient and optimum 
results have been obtained in typically eight 
to ten cycles. The results of the study 
indicate that there is an increase in the 
blade weight and a significant change in the 
design variable distributions with an increase 
in the number of frequency constraints. The 
optimization process tends to redistribute 
mass toward the blade tip due to the presence 
of the autorotational inertia constraint . The 
inclusion of the stress constraints has 
different effects on the wall thickness 
distributions of the rectangular and the 
tapered blades, but tends to increase the 
magnitude of the nonstructural segment weight 
distributions in both cases. 
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Table 1. Reference blade data (Fig. 2) 


Seg- 

Length 

Box beam dimension (ft) 

Bending stiff- 

Tor- 

Nonst ructural 

Pre- 

ment 

(ft) 


ness x 10 4 

s ional 

segment weight 

twist 

Num- 



(lb - ft 2 ) 

stiff- 

(lbs) 

(deg. ) 

ber 




ness 

w 

0 






4 

x 10* 







(Ib-ft 2 ) 




L 

t l t 2 C 3 

EI x EI y 

GJ 

Total Min. 



1 

1.37 

0.0116 

0.0080 

0.0280 

7.349 

78.58 

11 . Ill 

6.718 

0.89 

1 . 745 

2 

2.2 

0.0100 

0.0100 

0.0440 

6 . 957 

84.68 

10 .139 

9.088 

1 .435 

2.61’ 

3 

2.2 

0.0075 

0.0075 

0.0325 

5.548 

66.55 

7.778 

1.978 

1.435 

5.594 

4 

2.2 

0.0060 

0 . 0050 

0.0050 

4.128 

35 . 40 

5.833 

1.435 

1 . 435 

3.725 

5 

2.2 

0.0050 

0.0050 

0.0045 

3.537 

31.20 

5.000 

2 . 352 

1.435 

6 . 305 

6 

2.2 

0.0050 

0.0050 

0.0035 

3.514 

29.89 

4 .861 

5 .852 

1.435 

5.235 

7 

2.2 

0.0050 

0.0050 

0.0040 

3.526 

30.55 

4 .931 

6.342 

1.435 

3 .49 

8 

2.2 

0.0050 

0.0050 

0.0046 

3.539 

31 .31 

5.000 

6.573 

1.435 

0.00 

9 

2.2 

0.0050 

0 . 0050 

0 . 0035 

3.514 

29.39 

4.361 

6 . 372 

1 . 435 

-0.175 

10 

2.2 

0.0050 

0.0050 

0.0021 

3.481 

27 . 91 

2.778 

5.962 

1 . 435 

-1 . 915 


Table 2. Blade preassigned properties 


Number of blades 
Blade radius 


22 ft. 


Pitch axis droop 
Pitch axis sweep 
Rotor speed 


0 degree 
0 degree 
293 rpm 


Table 3. Reference blade frequencies and bounds 
(windows) 


Reference Blade 
Frequency 


Prescribed Bounds 
lower upper 


Chord 

1.3 ft. 


Hz 

per rev 

Hz 

per rev 

Hz 

per rev 

Flap hinge offset 

0.833 ft. 

f l 

12.295 

2.52 

12.162 

2.49 

12 . 408 

2.54 

Inplane hinge offset 

0.S33 ft. 

C 

“2 

16.098 

3.30 

15.936 

3.26 

16.258 

3.33 

Solidity (based on 
mean chord) 

0.0748 

e 

*■ 3 

20.913 

4.28 

20.704 

4.24 

21.122 

4.33 

Precone angle 

0 degree 

p 

*■ 4 

34 . 624 

7.09 

34.272 

7.02 

34 . 966 

7 .16 

Droop angle 
Tip sweep 

0 degree 
0 degree 

f 5 

35 .361 

7 .34 

35.502 

7.27 

36.219 

7 . 42 
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Table 4. Summary of cases studied 


Con- No. of Planform 

straint Design 
Case Variables 

Design 
variables 
(i-1,2, . . 10) 

a 

30 

Rectangular 

t l.' t 2.' t 3. 

Ill 

b 

30 

Rectangular 

t l. ,t 2. ,t 3. 

111 

c 

30 

Rectangular 

t. , ty f t ^ 
i i i 

a 

40 

Rectangular 

e l. ,t 2.' t 3.' W o. 
1111 

b 

40 

Rectangular 

t, ,t- ,t, ,w 

iiii 

c 

40 

Rectangular 

t l.' t 2. , '3. '% 
1111 

a 

42 

Tapered 

h , X. , t , , 1 0 , 
r' h' 1 . 2 . 

i i 

t, ,W 
3 i °i 

b 

42 

Tapered 

1 1 

t, , w 

Y °i 

c 

42 

Tapered 

h r' t l . ' t 2 . ' 




*3.'V 
1 1 

Case 

Constraint definition Abbreviation used 

a 

Windows on first three 
lead-lag and first two 
flapping frequencies, 
autorotaional inertia 
and stress constraints 

5 freq, AI,<J 

b 

Windows on first three 
lead-lag and first two 
flapping frequencies and 
autorotat ional inertia 
constraint 

5 freq, AI 

c 

windows on first lead-lag 
and first flapping freq. 
and autorotational inertia 
constraint (Ref. 19) 

2 freq,AI 


Table 5. Optimization results for rectangular 
blade; cases a-c, 30 design variables 
(see Table 4) 


Reference 

blade 


Opt imum 
blade 




5 Freq 
AI 
a 

case a 

5 Freq 
A I 

case b 

2 Freq 
AI 

case c 

f x (Hz) 

12.285 

12 . 408 

12 . 40 8 

12.438 

f 2 (Hz) 

15.098 

16 . 056 

16.344 

15.945 

f 3 (Hz) 

20 . 913 

20 . 968 

21.327 

20 . 37V 

f 4 (Hz' 

34.624 

34 . 546 

34.594 

33.363 

f 5 (Hz) 

35.361 

35 . 502 

35 . 502 

34.231 

Auto- 

rotational 

inertia 

( Ib-ft^) 

517 . 3 

517.3 

517.3 

517.3 

3lade weight 
(lb) 

98.27 

95 . 62 

95.23 

3 9.92 

Percent reduc 
tion in blade 

- 

2.67 

3.04 

3.50 


* 

weight 


* - From reference blade 


Table 6. Optimization results for rectangular 
blade; cases a-c, 40 design variables 


{ see 

Table 4) 




Reference 

blade 


Optimum 

blade 




5 Freq 
AI 
<J 

case a 

5 Freq 
AI 

case b 

2 Freq 
AI 

case c 

f x (Hz) 

12.285 

12 . 408 

12.408 

12.408 

f 2 (Hz) 

16.098 

16 . 075 

16.025 

15 . 940 

f 3 (Hz) 

20 .913 

21.081 

21.060 

22.600 

f 4 (Hz) 

34 . 624 

34 . 823 

34 . 589 

37 .350 

f 5 (Hz) 

35.361 

35.800 

35 . 595 

38.710 

Auto- 

rotational 

inertia 

517 . 3 

517 . 3 

517.3 

517.3 

(lb-ft 2 ) 





Blade weight 
(lb) 

98.27 

93.613 

90.624 

35.270 

Percent reduc- 
tion in blade 

• 

weight 

- 

4.74 

7.78 

13.23 


• - From reference blade 
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Table 7. Optimization results for tapered 

blade; cases a-c, 42 design variables 
(see Table 4) 


Reference 

blade 


Optimum 

blade 




5 Freq 
AI 

5 Freq 
AI 

2 Freq 
AI 



a 

- 

- 



case a 

case b 

case c 

X h 

1.0 

1 . 490 

1.508 

i . in 

f x (Hz) 

12.285 

12.408 

12.408 

12 . 408 

f 2 (Hz) 

16.098 

16.066 

16.064 

15.938 

f 3 (Hz) 

20 . 913 

20 . 888 

20.959 

22 . 504 

f 4 (Hz) 

34.624 

34 . 678 

34 . 646 

36.753 

f 5 (Hz) 

35.861 

35.507 

35 . 525 

38 .447 

AutO- 

rotational 

inertia 

517.3 

517 . 3 

517.3 

517 . 3 

(lb-ft 2 ) 





Blade weight 
(lb) 

98.27 

92 .16 

89.24 

84.24 

Percent reduc- 
tion in blade 

- 

6.21 

9.19 

14.28 

weight 






• - From reference blade 



Fig. 1 Flowchart of the optimization process 



Fig. 2 Rotor blade cross section 




z 


h(z) = h r (1 -z/R) + h, z/R 


b) Rotor blade taper 
Fig. 3 Rotor blade 
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008 

Horizontal 

wall thickness (t-f), .006 
ft 

.004 


Reference blade 
5 frequency, Al, a 



83 4.4 


8.8 13.2 17.6 22.0 

Blade radius, ft 


a) Rectangular blade, 30 design variables 

.012 r- Reference blade 

5 frequency, Al, a 

.010 b 1 1 


008 j- 

Horizonta! 

wall thickness (t-j ). .006 j- 
ft 1 

.004 b 


Nonstructural 
segment weight, 4.8 - 
lb 

3.2 - 



■ Reference blade 

■ 5 frequency. Al. o 



83 44 8.8 13.2 175 22.0 

Blade radius, ft 

a) Rectangular blade, 40 design variables 


Nonstructural 
segment weight. 4 8 
lb 



Reference blade 
5 frequency Al. q 


.83 4.4 8.8 13.2 1 7 6 22.0 

Blade radius, ft 

b) Tapered blade, 42 design variables 

Fig. 4 Optimum distribution of box beam 
horizontal wall thickness ^t . 1 ) along 
blade radius 


Reference blade 
5 frequency, Al, q 



Vertical 

wall thickness (( 2 ). 006 

ft 


.83 4.4 8.8 13.2 17.6 22.0 

Blade radius, ft 

a) Rectangular blade, 30 design variables 



Vertical 

wall thickness (t 2 >. .006 
ft 


.83 4.4 8.8 13.2 17.6 22.0 

Blade radius, ft 

b) Tapered blade, 42 design variables 

Fig. 5 Optimum distribution of box beam 
vertical wall thickness (tj) along 
blade radius 


.83 4.4 8.8 1 3.2 '7 5 220 

Blade radius, ft 

b) Tapered blade, 42 design variables 

Fig. 5 Optimum distribution of nonstructural 
segment weight along blade radius 

Reference blade 

Q12 5 frequency. Al 

j 2 frequency. Al 


Horizontal — i 

wall thickness (t^). .006 u | I 1 

ft ; --'l T= P 

004 r ! 1 Ur" — ' — 


83 4.4 8.8 1 3 2 17 6 22.0 

Blade radius, ft 

a) Rectangular blade, 30 design variables 

018 - 

Reference blade i ] 

01 5- 5 frequency. Al j j 

2 ‘requency. Al I 

012 — — * 

Horizontal | 1 ! 

wall thickness (tj). .009 ^ ^ I j j 

003 - l.j — I I 

o : 1 : : ! 


L.J 

83 4.4 


8.8 13.2 17 6 22.0 

Blade radius, ft 


b) Tapered blade, 42 design variables 

Fig. 7 Optimum distributions of box beam 
horizontal wall thickness (p) along 

blade radius; effect of higher 

frequency constraints 
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- a) Rectangular blade, 30 design variables a) Rectangular blade, 30 design variables 


Reference blade 


Vertical 
wall thickness 
ft 


b) 



Tapered blade, 42 design variables 


Fig. 3 Optimum distributions of box beam 
vertical wall thickness (t 2 ) along 

blade radius; effect of higher 

frequency constraints 


Reference blade 

5 frequency. A1 o 

5 frequency. A! 



Blade radius, ft 


b) Tapered blade, 42 design variables 

Fig. 10 Optimum distributions of box beam 
horizontal wall thickness ( t ^ ) along 

blade radius; effect of stress 
constraints 


Ncnstructural 
segment weight, 
lb 



Nonstructural 
segment weight, 
lb 



a) Rectangular blade, 40 design variables 


a) Rectangular blade, 40 design variables 



b) Tapered blade, 42 design variables 

ig. 9 Optimum distributions of nonstructural 
segment weight along blade radius 



b) Tapered blade, 42 design variables 

Fig. 11 Optimum distribution of nonstructural 
segment weight along blade radius; 
effect of stress constraints 
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